if (!require("dplyr")) install.packages("dplyr")
## Loading required package: dplyr
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
if (!require("skimr")) install.packages("skimr")
## Loading required package: skimr
if (!require("tidyr")) install.packages("tidyr")
## Loading required package: tidyr
if (!require("survival")) install.packages("survival")
## Loading required package: survival
if (!require("survminer")) install.packages("survminer")
## Loading required package: survminer
## Loading required package: ggplot2
## Loading required package: ggpubr
##
## Attaching package: 'survminer'
## The following object is masked from 'package:survival':
##
## myeloma
if (!require("haven")) install.packages("haven")
## Loading required package: haven
if (!require("broom")) install.packages("broom")
## Loading required package: broom
if (!require("rms")) install.packages("rms")
## Loading required package: rms
## Loading required package: Hmisc
##
## Attaching package: 'Hmisc'
## The following objects are masked from 'package:dplyr':
##
## src, summarize
## The following objects are masked from 'package:base':
##
## format.pval, units
if (!require("tidyverse")) install.packages("tidyverse")
## Loading required package: tidyverse
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats 1.0.0 ✔ readr 2.1.5
## ✔ lubridate 1.9.4 ✔ stringr 1.5.1
## ✔ purrr 1.0.4 ✔ tibble 3.2.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ✖ Hmisc::src() masks dplyr::src()
## ✖ Hmisc::summarize() masks dplyr::summarize()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
if (!require("tableone")) install.packages("tableone")
## Loading required package: tableone
library(dplyr)
library(skimr)
library(tidyr)
library(survival)
library(survminer)
library(haven)
library(broom)
library(rms)
library(tidyverse)
library(tableone)
NHANES2 <- read.csv("NHANES2-1 (1).csv")
d <- NHANES2 #%>%
#select('ROWNAMES','SEX','RACE','MARRY','DEATH','AGEYRS',
#'GRADES','WT', 'BOOZE', 'SIZE',
#'AVGSMK', "HEIGHT", "EXAM_YR", "DIE_YR", "LAST_YR")
#Exclude missing death
d <- d %>%
filter(!is.na(BOOZE), !is.na(DEATH), !is.na(SEX), !is.na(RACE), !is.na(GRADES), !is.na(MARRY), !is.na(AVGSMK), !is.na(SIZE), !is.na(GRADES))
#BMI
d <- d %>%
mutate(BMI = WT / (HEIGHT / 100)^2)
head(d$BMI)
## [1] 20.49522 21.02151 23.22748 35.72785 27.92312 30.50132
# GRADES and SIZE categories
d$EDUC_CAT <- cut(d$GRADES,
breaks = c(-Inf, 8, 11, 12, 15, Inf),
labels = c("≤8 yrs", "Some HS", "HS Grad", "Some College", "College+"),
right = TRUE)
d$SIZE_CAT <- cut(d$SIZE,
breaks = c(0, 3, 5, 7, 8),
labels = c("Rural", "Small town", "Medium city", "Large city"),
right = TRUE)
# Catergorical BOOZE
d <- d %>%
mutate(BOOZE_q = cut(
BOOZE,
breaks = c(-1, 0, 0.5, 2.0, 77.0),
include.lowest = TRUE,
labels = c("0/week", "0–0.5/week", "0.5–2/week", ">2/week")
))
vars <- c("AGEYRS", "SEX", "RACE", "MARRY", "BMI", "AVGSMK", "EDUC_CAT", "SIZE_CAT")
catVars <- c("SEX", "RACE", "MARRY")
#Table 1
table1 <- CreateTableOne(vars = vars,
data = d,
strata = "BOOZE_q",
factorVars = catVars)
print(table1, showAllLevels = TRUE)
## Stratified by BOOZE_q
## level 0/week 0–0.5/week 0.5–2/week
## n 4053 941 1729
## AGEYRS (mean (SD)) 57.09 (12.79) 54.34 (13.36) 51.60 (13.53)
## SEX (%) 1 1448 (35.7) 367 (39.0) 856 (49.5)
## 2 2605 (64.3) 574 (61.0) 873 (50.5)
## RACE (%) 1 3497 (86.3) 827 (87.9) 1515 (87.6)
## 2 475 (11.7) 93 ( 9.9) 194 (11.2)
## 3 81 ( 2.0) 21 ( 2.2) 20 ( 1.2)
## MARRY (%) 2 2885 (71.2) 683 (72.6) 1288 (74.5)
## 3 671 (16.6) 127 (13.5) 172 ( 9.9)
## 4 190 ( 4.7) 67 ( 7.1) 102 ( 5.9)
## 5 96 ( 2.4) 20 ( 2.1) 59 ( 3.4)
## 6 202 ( 5.0) 40 ( 4.3) 103 ( 6.0)
## 8 9 ( 0.2) 4 ( 0.4) 5 ( 0.3)
## BMI (mean (SD)) 26.55 (5.50) 26.42 (5.10) 26.03 (4.84)
## AVGSMK (mean (SD)) 4.82 (10.84) 6.72 (12.63) 8.26 (13.77)
## EDUC_CAT (%) ≤8 yrs 1452 (35.8) 217 (23.1) 337 (19.5)
## Some HS 753 (18.6) 185 (19.7) 282 (16.3)
## HS Grad 1209 (29.8) 311 (33.0) 636 (36.8)
## Some College 353 ( 8.7) 130 (13.8) 235 (13.6)
## College+ 286 ( 7.1) 98 (10.4) 239 (13.8)
## SIZE_CAT (%) Rural 1101 (27.2) 348 (37.0) 758 (43.8)
## Small town 454 (11.2) 123 (13.1) 244 (14.1)
## Medium city 569 (14.0) 118 (12.5) 205 (11.9)
## Large city 1929 (47.6) 352 (37.4) 522 (30.2)
## Stratified by BOOZE_q
## >2/week p test
## n 2527
## AGEYRS (mean (SD)) 51.71 (13.18) <0.001
## SEX (%) 1678 (66.4) <0.001
## 849 (33.6)
## RACE (%) 2250 (89.0) 0.010
## 236 ( 9.3)
## 41 ( 1.6)
## MARRY (%) 1972 (78.0) <0.001
## 170 ( 6.7)
## 168 ( 6.6)
## 70 ( 2.8)
## 140 ( 5.5)
## 7 ( 0.3)
## BMI (mean (SD)) 25.20 (4.08) <0.001
## AVGSMK (mean (SD)) 9.60 (13.92) <0.001
## EDUC_CAT (%) 387 (15.3) <0.001
## 354 (14.0)
## 876 (34.7)
## 411 (16.3)
## 499 (19.7)
## SIZE_CAT (%) 1239 (49.0) <0.001
## 346 (13.7)
## 261 (10.3)
## 681 (26.9)
#Create follow-up time
d$start <- d$EXAM_YR + d$EXAM_MO / 12
d$end <- ifelse(d$DEATH == 1,
d$DIE_YR + d$DIE_MO / 12,
d$LAST_YR + d$LAST_MO / 12)
d$FU <- d$end - d$start
#Adjusted Cox regression with SEX
cox <- coxph(Surv(FU, DEATH) ~ as.factor(BOOZE_q) + SEX +
as.factor(RACE) + as.factor(GRADES) + as.factor(MARRY) + BMI +
AVGSMK + SIZE, data = d, ties='efron')
summary(cox)
## Call:
## coxph(formula = Surv(FU, DEATH) ~ as.factor(BOOZE_q) + SEX +
## as.factor(RACE) + as.factor(GRADES) + as.factor(MARRY) +
## BMI + AVGSMK + SIZE, data = d, ties = "efron")
##
## n= 9250, number of events= 2145
##
## coef exp(coef) se(coef) z Pr(>|z|)
## as.factor(BOOZE_q)0–0.5/week -0.131786 0.876529 0.075867 -1.737 0.082376 .
## as.factor(BOOZE_q)0.5–2/week -0.390542 0.676690 0.065568 -5.956 2.58e-09 ***
## as.factor(BOOZE_q)>2/week -0.333146 0.716665 0.058474 -5.697 1.22e-08 ***
## SEX -0.773141 0.461561 0.049172 -15.723 < 2e-16 ***
## as.factor(RACE)2 -0.265985 0.766450 0.075590 -3.519 0.000434 ***
## as.factor(RACE)3 -0.513794 0.598222 0.198879 -2.583 0.009782 **
## as.factor(GRADES)1 0.599987 1.822095 0.314018 1.911 0.056046 .
## as.factor(GRADES)2 0.520093 1.682184 0.295838 1.758 0.078742 .
## as.factor(GRADES)3 0.701697 2.017172 0.252651 2.777 0.005481 **
## as.factor(GRADES)4 0.668808 1.951909 0.244717 2.733 0.006276 **
## as.factor(GRADES)5 0.707530 2.028973 0.250421 2.825 0.004723 **
## as.factor(GRADES)6 0.344214 1.410881 0.236150 1.458 0.144949
## as.factor(GRADES)7 0.391113 1.478625 0.237484 1.647 0.099578 .
## as.factor(GRADES)8 0.524292 1.689263 0.220388 2.379 0.017362 *
## as.factor(GRADES)9 0.393780 1.482575 0.230085 1.711 0.086997 .
## as.factor(GRADES)10 0.177886 1.194689 0.228746 0.778 0.436771
## as.factor(GRADES)11 0.214700 1.239490 0.235460 0.912 0.361857
## as.factor(GRADES)12 -0.033569 0.966988 0.218656 -0.154 0.877984
## as.factor(GRADES)13 -0.132898 0.875555 0.243908 -0.545 0.585843
## as.factor(GRADES)14 -0.214628 0.806841 0.241954 -0.887 0.375046
## as.factor(GRADES)15 -0.177631 0.837251 0.272086 -0.653 0.513853
## as.factor(GRADES)16 -0.170102 0.843579 0.237427 -0.716 0.473721
## as.factor(GRADES)17 -0.615065 0.540606 0.248169 -2.478 0.013197 *
## as.factor(MARRY)3 0.708942 2.031840 0.062116 11.413 < 2e-16 ***
## as.factor(MARRY)4 0.065615 1.067815 0.101914 0.644 0.519686
## as.factor(MARRY)5 0.030171 1.030631 0.144768 0.208 0.834910
## as.factor(MARRY)6 0.156516 1.169430 0.096475 1.622 0.104728
## as.factor(MARRY)8 0.746167 2.108901 0.336343 2.218 0.026523 *
## BMI -0.013787 0.986307 0.004682 -2.945 0.003232 **
## AVGSMK 0.004971 1.004984 0.001629 3.051 0.002279 **
## SIZE -0.022890 0.977370 0.008506 -2.691 0.007120 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## as.factor(BOOZE_q)0–0.5/week 0.8765 1.1409 0.7554 1.0171
## as.factor(BOOZE_q)0.5–2/week 0.6767 1.4778 0.5951 0.7695
## as.factor(BOOZE_q)>2/week 0.7167 1.3954 0.6391 0.8037
## SEX 0.4616 2.1666 0.4192 0.5083
## as.factor(RACE)2 0.7665 1.3047 0.6609 0.8888
## as.factor(RACE)3 0.5982 1.6716 0.4051 0.8834
## as.factor(GRADES)1 1.8221 0.5488 0.9846 3.3718
## as.factor(GRADES)2 1.6822 0.5945 0.9420 3.0039
## as.factor(GRADES)3 2.0172 0.4957 1.2294 3.3098
## as.factor(GRADES)4 1.9519 0.5123 1.2082 3.1533
## as.factor(GRADES)5 2.0290 0.4929 1.2420 3.3146
## as.factor(GRADES)6 1.4109 0.7088 0.8881 2.2413
## as.factor(GRADES)7 1.4786 0.6763 0.9283 2.3551
## as.factor(GRADES)8 1.6893 0.5920 1.0967 2.6019
## as.factor(GRADES)9 1.4826 0.6745 0.9444 2.3274
## as.factor(GRADES)10 1.1947 0.8370 0.7630 1.8705
## as.factor(GRADES)11 1.2395 0.8068 0.7813 1.9664
## as.factor(GRADES)12 0.9670 1.0341 0.6299 1.4844
## as.factor(GRADES)13 0.8756 1.1421 0.5428 1.4122
## as.factor(GRADES)14 0.8068 1.2394 0.5022 1.2964
## as.factor(GRADES)15 0.8373 1.1944 0.4912 1.4271
## as.factor(GRADES)16 0.8436 1.1854 0.5297 1.3435
## as.factor(GRADES)17 0.5406 1.8498 0.3324 0.8793
## as.factor(MARRY)3 2.0318 0.4922 1.7989 2.2949
## as.factor(MARRY)4 1.0678 0.9365 0.8745 1.3039
## as.factor(MARRY)5 1.0306 0.9703 0.7760 1.3688
## as.factor(MARRY)6 1.1694 0.8551 0.9680 1.4128
## as.factor(MARRY)8 2.1089 0.4742 1.0908 4.0771
## BMI 0.9863 1.0139 0.9773 0.9954
## AVGSMK 1.0050 0.9950 1.0018 1.0082
## SIZE 0.9774 1.0232 0.9612 0.9938
##
## Concordance= 0.661 (se = 0.006 )
## Likelihood ratio test= 651 on 31 df, p=<2e-16
## Wald test = 649.1 on 31 df, p=<2e-16
## Score (logrank) test = 673.9 on 31 df, p=<2e-16
cox.zph(cox)
## chisq df p
## as.factor(BOOZE_q) 7.813 3 0.05004
## SEX 11.185 1 0.00082
## as.factor(RACE) 1.933 2 0.38039
## as.factor(GRADES) 20.679 17 0.24099
## as.factor(MARRY) 4.452 5 0.48634
## BMI 0.173 1 0.67753
## AVGSMK 0.458 1 0.49844
## SIZE 0.264 1 0.60707
## GLOBAL 48.114 31 0.02566
plot(cox.zph(cox))








#Product term with SEX Cox regression
cox_product <- coxph(Surv(FU, DEATH) ~ as.factor(BOOZE_q)*SEX +
as.factor(RACE) + as.factor(GRADES) + as.factor(MARRY) +
BMI + AVGSMK + SIZE, data = d, ties = "efron")
summary(cox_product)
## Call:
## coxph(formula = Surv(FU, DEATH) ~ as.factor(BOOZE_q) * SEX +
## as.factor(RACE) + as.factor(GRADES) + as.factor(MARRY) +
## BMI + AVGSMK + SIZE, data = d, ties = "efron")
##
## n= 9250, number of events= 2145
##
## coef exp(coef) se(coef) z Pr(>|z|)
## as.factor(BOOZE_q)0–0.5/week -0.161405 0.850947 0.236705 -0.682 0.495313
## as.factor(BOOZE_q)0.5–2/week -0.453902 0.635145 0.194085 -2.339 0.019352
## as.factor(BOOZE_q)>2/week -0.567586 0.566892 0.166506 -3.409 0.000653
## SEX -0.817512 0.441529 0.064487 -12.677 < 2e-16
## as.factor(RACE)2 -0.266209 0.766279 0.075599 -3.521 0.000429
## as.factor(RACE)3 -0.511089 0.599842 0.198899 -2.570 0.010182
## as.factor(GRADES)1 0.597777 1.818073 0.314051 1.903 0.056983
## as.factor(GRADES)2 0.516572 1.676272 0.295965 1.745 0.080919
## as.factor(GRADES)3 0.705668 2.025200 0.252666 2.793 0.005224
## as.factor(GRADES)4 0.674080 1.962227 0.244765 2.754 0.005887
## as.factor(GRADES)5 0.712374 2.038825 0.250470 2.844 0.004453
## as.factor(GRADES)6 0.347922 1.416121 0.236162 1.473 0.140688
## as.factor(GRADES)7 0.394282 1.483318 0.237521 1.660 0.096917
## as.factor(GRADES)8 0.527903 1.695374 0.220403 2.395 0.016612
## as.factor(GRADES)9 0.394902 1.484239 0.230093 1.716 0.086112
## as.factor(GRADES)10 0.181665 1.199212 0.228753 0.794 0.427108
## as.factor(GRADES)11 0.217739 1.243262 0.235479 0.925 0.355141
## as.factor(GRADES)12 -0.032482 0.968040 0.218671 -0.149 0.881915
## as.factor(GRADES)13 -0.134216 0.874401 0.243941 -0.550 0.582183
## as.factor(GRADES)14 -0.216948 0.804972 0.241980 -0.897 0.369958
## as.factor(GRADES)15 -0.178067 0.836886 0.272107 -0.654 0.512854
## as.factor(GRADES)16 -0.168185 0.845198 0.237443 -0.708 0.478749
## as.factor(GRADES)17 -0.613072 0.541684 0.248185 -2.470 0.013503
## as.factor(MARRY)3 0.713826 2.041788 0.062203 11.476 < 2e-16
## as.factor(MARRY)4 0.063219 1.065260 0.101959 0.620 0.535231
## as.factor(MARRY)5 0.032512 1.033046 0.144813 0.225 0.822363
## as.factor(MARRY)6 0.158054 1.171229 0.096484 1.638 0.101392
## as.factor(MARRY)8 0.747464 2.111638 0.336519 2.221 0.026340
## BMI -0.013141 0.986945 0.004706 -2.793 0.005227
## AVGSMK 0.004964 1.004976 0.001629 3.046 0.002316
## SIZE -0.022786 0.977472 0.008512 -2.677 0.007428
## as.factor(BOOZE_q)0–0.5/week:SEX 0.019546 1.019738 0.150620 0.130 0.896748
## as.factor(BOOZE_q)0.5–2/week:SEX 0.042411 1.043323 0.131314 0.323 0.746717
## as.factor(BOOZE_q)>2/week:SEX 0.183497 1.201412 0.120545 1.522 0.127949
##
## as.factor(BOOZE_q)0–0.5/week
## as.factor(BOOZE_q)0.5–2/week *
## as.factor(BOOZE_q)>2/week ***
## SEX ***
## as.factor(RACE)2 ***
## as.factor(RACE)3 *
## as.factor(GRADES)1 .
## as.factor(GRADES)2 .
## as.factor(GRADES)3 **
## as.factor(GRADES)4 **
## as.factor(GRADES)5 **
## as.factor(GRADES)6
## as.factor(GRADES)7 .
## as.factor(GRADES)8 *
## as.factor(GRADES)9 .
## as.factor(GRADES)10
## as.factor(GRADES)11
## as.factor(GRADES)12
## as.factor(GRADES)13
## as.factor(GRADES)14
## as.factor(GRADES)15
## as.factor(GRADES)16
## as.factor(GRADES)17 *
## as.factor(MARRY)3 ***
## as.factor(MARRY)4
## as.factor(MARRY)5
## as.factor(MARRY)6
## as.factor(MARRY)8 *
## BMI **
## AVGSMK **
## SIZE **
## as.factor(BOOZE_q)0–0.5/week:SEX
## as.factor(BOOZE_q)0.5–2/week:SEX
## as.factor(BOOZE_q)>2/week:SEX
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## as.factor(BOOZE_q)0–0.5/week 0.8509 1.1752 0.5351 1.3533
## as.factor(BOOZE_q)0.5–2/week 0.6351 1.5744 0.4342 0.9291
## as.factor(BOOZE_q)>2/week 0.5669 1.7640 0.4090 0.7857
## SEX 0.4415 2.2649 0.3891 0.5010
## as.factor(RACE)2 0.7663 1.3050 0.6607 0.8887
## as.factor(RACE)3 0.5998 1.6671 0.4062 0.8858
## as.factor(GRADES)1 1.8181 0.5500 0.9824 3.3646
## as.factor(GRADES)2 1.6763 0.5966 0.9385 2.9941
## as.factor(GRADES)3 2.0252 0.4938 1.2342 3.3231
## as.factor(GRADES)4 1.9622 0.5096 1.2145 3.1703
## as.factor(GRADES)5 2.0388 0.4905 1.2479 3.3310
## as.factor(GRADES)6 1.4161 0.7062 0.8914 2.2497
## as.factor(GRADES)7 1.4833 0.6742 0.9312 2.3627
## as.factor(GRADES)8 1.6954 0.5898 1.1007 2.6114
## as.factor(GRADES)9 1.4842 0.6737 0.9455 2.3300
## as.factor(GRADES)10 1.1992 0.8339 0.7659 1.8776
## as.factor(GRADES)11 1.2433 0.8043 0.7837 1.9724
## as.factor(GRADES)12 0.9680 1.0330 0.6306 1.4860
## as.factor(GRADES)13 0.8744 1.1436 0.5421 1.4104
## as.factor(GRADES)14 0.8050 1.2423 0.5010 1.2935
## as.factor(GRADES)15 0.8369 1.1949 0.4910 1.4265
## as.factor(GRADES)16 0.8452 1.1832 0.5307 1.3461
## as.factor(GRADES)17 0.5417 1.8461 0.3330 0.8811
## as.factor(MARRY)3 2.0418 0.4898 1.8074 2.3065
## as.factor(MARRY)4 1.0653 0.9387 0.8723 1.3009
## as.factor(MARRY)5 1.0330 0.9680 0.7778 1.3721
## as.factor(MARRY)6 1.1712 0.8538 0.9694 1.4150
## as.factor(MARRY)8 2.1116 0.4736 1.0919 4.0838
## BMI 0.9869 1.0132 0.9779 0.9961
## AVGSMK 1.0050 0.9950 1.0018 1.0082
## SIZE 0.9775 1.0230 0.9613 0.9939
## as.factor(BOOZE_q)0–0.5/week:SEX 1.0197 0.9806 0.7591 1.3699
## as.factor(BOOZE_q)0.5–2/week:SEX 1.0433 0.9585 0.8066 1.3496
## as.factor(BOOZE_q)>2/week:SEX 1.2014 0.8324 0.9486 1.5216
##
## Concordance= 0.661 (se = 0.006 )
## Likelihood ratio test= 653.3 on 34 df, p=<2e-16
## Wald test = 658.5 on 34 df, p=<2e-16
## Score (logrank) test = 694.2 on 34 df, p=<2e-16
cox.zph(cox_product)
## chisq df p
## as.factor(BOOZE_q) 7.669 3 0.05337
## SEX 11.045 1 0.00089
## as.factor(RACE) 1.948 2 0.37762
## as.factor(GRADES) 20.607 17 0.24437
## as.factor(MARRY) 4.475 5 0.48325
## BMI 0.172 1 0.67826
## AVGSMK 0.472 1 0.49207
## SIZE 0.253 1 0.61527
## as.factor(BOOZE_q):SEX 12.584 3 0.00563
## GLOBAL 47.959 34 0.05671
plot(cox.zph(cox_product))









#Kaplan
fit<-survfit(Surv(FU, DEATH)~BOOZE_q, data=d)
summary(fit)
## Call: survfit(formula = Surv(FU, DEATH) ~ BOOZE_q, data = d)
##
## BOOZE_q=0/week
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0000 4053 2 1.000 0.000349 0.999 1.000
## 0.0833 4051 3 0.999 0.000551 0.998 1.000
## 0.1667 4048 1 0.999 0.000604 0.997 1.000
## 0.2500 4047 4 0.998 0.000779 0.996 0.999
## 0.3333 4043 6 0.996 0.000985 0.994 0.998
## 0.4167 4037 2 0.996 0.001044 0.994 0.998
## 0.5000 4035 2 0.995 0.001101 0.993 0.997
## 0.5833 4033 7 0.993 0.001278 0.991 0.996
## 0.6667 4026 3 0.993 0.001346 0.990 0.995
## 0.7500 4023 3 0.992 0.001412 0.989 0.995
## 0.8333 4020 4 0.991 0.001494 0.988 0.994
## 0.9167 4016 1 0.991 0.001514 0.988 0.994
## 1.0000 4015 6 0.989 0.001628 0.986 0.992
## 1.0833 4009 5 0.988 0.001717 0.985 0.991
## 1.1667 4004 5 0.987 0.001801 0.983 0.990
## 1.2500 3999 5 0.985 0.001881 0.982 0.989
## 1.3333 3994 4 0.984 0.001943 0.981 0.988
## 1.4167 3990 4 0.983 0.002003 0.980 0.987
## 1.5000 3986 6 0.982 0.002089 0.978 0.986
## 1.5833 3980 4 0.981 0.002144 0.977 0.985
## 1.6667 3976 6 0.980 0.002225 0.975 0.984
## 1.7500 3970 10 0.977 0.002352 0.972 0.982
## 1.8333 3960 7 0.975 0.002437 0.971 0.980
## 1.9167 3953 6 0.974 0.002507 0.969 0.979
## 2.0000 3947 4 0.973 0.002552 0.968 0.978
## 2.0833 3943 4 0.972 0.002597 0.967 0.977
## 2.1667 3939 1 0.972 0.002608 0.967 0.977
## 2.2500 3938 7 0.970 0.002684 0.965 0.975
## 2.3333 3931 3 0.969 0.002716 0.964 0.974
## 2.4167 3928 7 0.967 0.002788 0.962 0.973
## 2.5000 3921 11 0.965 0.002898 0.959 0.970
## 2.5833 3910 8 0.963 0.002975 0.957 0.969
## 2.6667 3902 8 0.961 0.003050 0.955 0.967
## 2.7500 3894 2 0.960 0.003068 0.954 0.966
## 2.8333 3892 2 0.960 0.003086 0.954 0.966
## 2.9167 3890 7 0.958 0.003149 0.952 0.964
## 3.0000 3883 10 0.956 0.003236 0.949 0.962
## 3.0833 3873 6 0.954 0.003287 0.948 0.961
## 3.1667 3867 5 0.953 0.003329 0.946 0.959
## 3.2500 3862 6 0.951 0.003378 0.945 0.958
## 3.3333 3856 1 0.951 0.003386 0.945 0.958
## 3.4167 3855 6 0.950 0.003434 0.943 0.956
## 3.5000 3849 1 0.949 0.003442 0.943 0.956
## 3.5833 3848 5 0.948 0.003482 0.941 0.955
## 3.6667 3843 8 0.946 0.003544 0.939 0.953
## 3.7500 3835 10 0.944 0.003619 0.937 0.951
## 3.8333 3825 1 0.943 0.003627 0.936 0.951
## 3.9167 3824 6 0.942 0.003671 0.935 0.949
## 4.0000 3818 10 0.940 0.003743 0.932 0.947
## 4.0833 3808 8 0.938 0.003800 0.930 0.945
## 4.1667 3800 1 0.937 0.003807 0.930 0.945
## 4.2500 3799 3 0.937 0.003828 0.929 0.944
## 4.3333 3796 5 0.935 0.003862 0.928 0.943
## 4.4167 3791 4 0.934 0.003890 0.927 0.942
## 4.5000 3787 4 0.933 0.003917 0.926 0.941
## 4.5833 3783 5 0.932 0.003950 0.924 0.940
## 4.6667 3778 7 0.930 0.003997 0.923 0.938
## 4.7500 3771 5 0.929 0.004029 0.921 0.937
## 4.8333 3766 7 0.927 0.004074 0.920 0.935
## 4.9167 3759 9 0.925 0.004131 0.917 0.933
## 5.0000 3750 4 0.924 0.004156 0.916 0.932
## 5.0833 3746 6 0.923 0.004193 0.915 0.931
## 5.1667 3740 7 0.921 0.004236 0.913 0.929
## 5.2500 3733 5 0.920 0.004266 0.911 0.928
## 5.3333 3728 5 0.919 0.004296 0.910 0.927
## 5.4167 3723 7 0.917 0.004337 0.908 0.925
## 5.5000 3716 8 0.915 0.004383 0.906 0.924
## 5.5833 3708 6 0.913 0.004418 0.905 0.922
## 5.6667 3702 9 0.911 0.004469 0.902 0.920
## 5.7500 3693 9 0.909 0.004519 0.900 0.918
## 5.8333 3684 9 0.907 0.004568 0.898 0.916
## 5.9167 3675 8 0.905 0.004611 0.896 0.914
## 6.0000 3667 10 0.902 0.004664 0.893 0.911
## 6.0833 3657 6 0.901 0.004695 0.892 0.910
## 6.1667 3651 4 0.900 0.004716 0.891 0.909
## 6.2500 3647 4 0.899 0.004736 0.890 0.908
## 6.3333 3643 6 0.897 0.004767 0.888 0.907
## 6.4167 3637 6 0.896 0.004797 0.887 0.905
## 6.5000 3631 4 0.895 0.004817 0.886 0.904
## 6.5833 3627 6 0.893 0.004847 0.884 0.903
## 6.6667 3621 10 0.891 0.004896 0.881 0.901
## 6.7500 3611 5 0.890 0.004920 0.880 0.899
## 6.8333 3606 9 0.887 0.004963 0.878 0.897
## 6.9167 3597 7 0.886 0.004997 0.876 0.896
## 7.0000 3590 4 0.885 0.005015 0.875 0.895
## 7.0833 3586 7 0.883 0.005048 0.873 0.893
## 7.1667 3579 7 0.881 0.005080 0.871 0.891
## 7.2500 3572 3 0.881 0.005094 0.871 0.891
## 7.3333 3569 7 0.879 0.005125 0.869 0.889
## 7.4167 3562 1 0.879 0.005130 0.869 0.889
## 7.5000 3561 11 0.876 0.005179 0.866 0.886
## 7.5833 3550 8 0.874 0.005214 0.864 0.884
## 7.6667 3542 9 0.872 0.005253 0.861 0.882
## 7.7500 3533 7 0.870 0.005283 0.860 0.880
## 7.8333 3526 6 0.868 0.005308 0.858 0.879
## 7.9167 3520 10 0.866 0.005350 0.856 0.877
## 8.0000 3510 7 0.864 0.005379 0.854 0.875
## 8.0833 3503 3 0.864 0.005392 0.853 0.874
## 8.1667 3500 8 0.862 0.005424 0.851 0.872
## 8.2500 3492 5 0.860 0.005445 0.850 0.871
## 8.3333 3487 5 0.859 0.005465 0.848 0.870
## 8.4167 3482 9 0.857 0.005500 0.846 0.868
## 8.5000 3473 7 0.855 0.005528 0.844 0.866
## 8.5833 3466 9 0.853 0.005563 0.842 0.864
## 8.6667 3457 6 0.851 0.005586 0.841 0.862
## 8.7500 3451 4 0.850 0.005601 0.840 0.862
## 8.8333 3447 4 0.849 0.005617 0.839 0.861
## 8.9167 3443 8 0.848 0.005647 0.837 0.859
## 9.0000 3435 6 0.846 0.005669 0.835 0.857
## 9.0833 3429 3 0.845 0.005680 0.834 0.857
## 9.1667 3426 10 0.843 0.005717 0.832 0.854
## 9.2500 3416 2 0.842 0.005724 0.831 0.854
## 9.3333 3414 4 0.841 0.005739 0.830 0.853
## 9.4167 3410 4 0.840 0.005753 0.829 0.852
## 9.5000 3406 4 0.839 0.005768 0.828 0.851
## 9.5833 3402 8 0.837 0.005796 0.826 0.849
## 9.6667 3394 8 0.835 0.005824 0.824 0.847
## 9.7500 3386 4 0.834 0.005838 0.823 0.846
## 9.8333 3382 10 0.832 0.005873 0.821 0.844
## 9.9167 3372 5 0.831 0.005890 0.819 0.842
## 10.0000 3367 8 0.829 0.005917 0.817 0.840
## 10.0833 3359 8 0.827 0.005944 0.815 0.839
## 10.1667 3351 11 0.824 0.005981 0.812 0.836
## 10.2500 3340 3 0.823 0.005991 0.812 0.835
## 10.3333 3337 10 0.821 0.006023 0.809 0.833
## 10.4167 3327 8 0.819 0.006049 0.807 0.831
## 10.5000 3319 6 0.817 0.006068 0.806 0.829
## 10.5833 3313 4 0.816 0.006081 0.805 0.828
## 10.6667 3309 7 0.815 0.006103 0.803 0.827
## 10.7500 3302 8 0.813 0.006128 0.801 0.825
## 10.8333 3294 5 0.811 0.006143 0.800 0.824
## 10.9167 3289 8 0.810 0.006168 0.798 0.822
## 11.0000 3281 6 0.808 0.006186 0.796 0.820
## 11.0833 3275 5 0.807 0.006201 0.795 0.819
## 11.1667 3270 6 0.805 0.006219 0.793 0.818
## 11.2500 3264 5 0.804 0.006234 0.792 0.816
## 11.3333 3259 12 0.801 0.006270 0.789 0.814
## 11.4167 3247 5 0.800 0.006284 0.788 0.812
## 11.5000 3242 10 0.797 0.006313 0.785 0.810
## 11.5833 3232 9 0.795 0.006339 0.783 0.808
## 11.6667 3223 9 0.793 0.006364 0.781 0.806
## 11.7500 3214 5 0.792 0.006378 0.779 0.804
## 11.8333 3209 5 0.791 0.006392 0.778 0.803
## 11.9167 3204 7 0.789 0.006411 0.776 0.801
## 12.0000 3197 6 0.787 0.006428 0.775 0.800
## 12.0833 3191 8 0.785 0.006449 0.773 0.798
## 12.1667 3183 7 0.784 0.006468 0.771 0.796
## 12.2500 3176 3 0.783 0.006476 0.770 0.796
## 12.3333 3173 9 0.781 0.006500 0.768 0.794
## 12.4167 3164 7 0.779 0.006518 0.766 0.792
## 12.5000 3157 9 0.777 0.006541 0.764 0.790
## 12.5833 3148 9 0.774 0.006565 0.762 0.787
## 12.6667 3139 3 0.774 0.006572 0.761 0.787
## 12.7500 3136 8 0.772 0.006592 0.759 0.785
## 12.8333 3128 6 0.770 0.006607 0.757 0.783
## 12.9167 3032 6 0.769 0.006624 0.756 0.782
## 13.0000 2965 9 0.766 0.006649 0.754 0.780
## 13.0833 2935 5 0.765 0.006663 0.752 0.778
## 13.1667 2849 4 0.764 0.006675 0.751 0.777
## 13.2500 2773 4 0.763 0.006689 0.750 0.776
## 13.3333 2688 6 0.761 0.006710 0.748 0.775
## 13.4167 2583 4 0.760 0.006725 0.747 0.773
## 13.5000 2578 6 0.758 0.006748 0.745 0.772
## 13.5833 2503 3 0.757 0.006760 0.744 0.771
## 13.6667 2438 7 0.755 0.006791 0.742 0.769
## 13.7500 2342 5 0.754 0.006814 0.740 0.767
## 13.8333 2194 5 0.752 0.006842 0.739 0.765
## 13.9167 2115 3 0.751 0.006860 0.737 0.764
## 14.0000 2062 1 0.750 0.006866 0.737 0.764
## 14.0833 2048 2 0.750 0.006879 0.736 0.763
## 14.1667 2009 1 0.749 0.006886 0.736 0.763
## 14.2500 1965 6 0.747 0.006928 0.734 0.761
## 14.3333 1926 6 0.745 0.006971 0.731 0.759
## 14.4167 1857 2 0.744 0.006987 0.730 0.758
## 14.5000 1806 6 0.741 0.007036 0.728 0.755
## 14.5833 1700 2 0.741 0.007055 0.727 0.755
## 14.6667 1617 4 0.739 0.007096 0.725 0.753
## 14.7500 1518 2 0.738 0.007120 0.724 0.752
## 14.8333 1456 2 0.737 0.007147 0.723 0.751
## 14.9167 1425 1 0.736 0.007160 0.722 0.750
## 15.0000 1397 4 0.734 0.007217 0.720 0.748
## 15.0833 1340 3 0.733 0.007263 0.718 0.747
## 15.1667 1291 2 0.731 0.007296 0.717 0.746
## 15.2500 1262 1 0.731 0.007313 0.717 0.745
## 15.3333 1228 5 0.728 0.007403 0.713 0.742
## 15.4167 1178 3 0.726 0.007461 0.711 0.741
## 15.5000 1145 2 0.725 0.007502 0.710 0.740
## 15.5833 1114 2 0.723 0.007545 0.709 0.738
## 15.6667 1079 2 0.722 0.007590 0.707 0.737
## 15.7500 1018 1 0.721 0.007616 0.707 0.736
## 15.8333 860 1 0.721 0.007653 0.706 0.736
## 15.9167 803 3 0.718 0.007781 0.703 0.733
## 16.0000 770 2 0.716 0.007871 0.701 0.732
## 16.0833 670 1 0.715 0.007932 0.699 0.731
## 16.2500 548 1 0.714 0.008024 0.698 0.729
## 16.3333 512 3 0.709 0.008332 0.693 0.726
## 16.5833 220 1 0.706 0.008896 0.689 0.724
##
## BOOZE_q=0–0.5/week
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0000 941 1 0.999 0.00106 0.997 1.000
## 0.0833 940 2 0.997 0.00184 0.993 1.000
## 0.2500 938 1 0.996 0.00212 0.992 1.000
## 0.3333 937 1 0.995 0.00237 0.990 0.999
## 0.5000 936 1 0.994 0.00259 0.989 0.999
## 0.5833 935 2 0.991 0.00299 0.986 0.997
## 0.6667 933 1 0.990 0.00317 0.984 0.997
## 0.8333 932 1 0.989 0.00334 0.983 0.996
## 0.9167 931 1 0.988 0.00350 0.981 0.995
## 1.1667 930 2 0.986 0.00381 0.979 0.994
## 1.2500 928 1 0.985 0.00395 0.977 0.993
## 1.3333 927 1 0.984 0.00408 0.976 0.992
## 1.4167 926 1 0.983 0.00421 0.975 0.991
## 1.5000 925 2 0.981 0.00447 0.972 0.990
## 1.7500 923 2 0.979 0.00470 0.970 0.988
## 1.8333 921 1 0.978 0.00482 0.968 0.987
## 2.0000 920 1 0.977 0.00493 0.967 0.986
## 2.0833 919 1 0.976 0.00503 0.966 0.985
## 2.1667 918 1 0.974 0.00514 0.964 0.985
## 2.2500 917 2 0.972 0.00534 0.962 0.983
## 2.3333 915 1 0.971 0.00544 0.961 0.982
## 2.4167 914 1 0.970 0.00554 0.959 0.981
## 2.5000 913 1 0.969 0.00563 0.958 0.980
## 2.5833 912 1 0.968 0.00573 0.957 0.979
## 2.6667 911 2 0.966 0.00591 0.954 0.978
## 2.7500 909 2 0.964 0.00608 0.952 0.976
## 2.8333 907 2 0.962 0.00625 0.950 0.974
## 2.9167 905 1 0.961 0.00634 0.948 0.973
## 3.0000 904 1 0.960 0.00642 0.947 0.972
## 3.0833 903 1 0.959 0.00650 0.946 0.971
## 3.2500 902 3 0.955 0.00673 0.942 0.969
## 3.3333 899 2 0.953 0.00688 0.940 0.967
## 3.5000 897 1 0.952 0.00696 0.939 0.966
## 3.6667 896 2 0.950 0.00710 0.936 0.964
## 3.7500 894 3 0.947 0.00731 0.933 0.961
## 3.8333 891 1 0.946 0.00738 0.931 0.960
## 4.0000 890 1 0.945 0.00745 0.930 0.959
## 4.0833 889 1 0.944 0.00752 0.929 0.959
## 4.4167 888 1 0.943 0.00758 0.928 0.958
## 4.5000 887 1 0.942 0.00765 0.927 0.957
## 4.7500 886 3 0.938 0.00784 0.923 0.954
## 4.9167 883 1 0.937 0.00790 0.922 0.953
## 5.0000 882 1 0.936 0.00796 0.921 0.952
## 5.0833 881 2 0.934 0.00809 0.918 0.950
## 5.3333 879 2 0.932 0.00821 0.916 0.948
## 5.4167 877 1 0.931 0.00827 0.915 0.947
## 5.5833 876 1 0.930 0.00833 0.914 0.946
## 5.6667 875 1 0.929 0.00838 0.913 0.945
## 5.7500 874 2 0.927 0.00850 0.910 0.943
## 5.8333 872 1 0.926 0.00855 0.909 0.943
## 5.9167 871 1 0.925 0.00861 0.908 0.942
## 6.1667 870 4 0.920 0.00883 0.903 0.938
## 6.2500 866 1 0.919 0.00888 0.902 0.937
## 6.3333 865 2 0.917 0.00899 0.900 0.935
## 6.5000 863 1 0.916 0.00904 0.898 0.934
## 6.5833 862 3 0.913 0.00919 0.895 0.931
## 6.6667 859 1 0.912 0.00924 0.894 0.930
## 7.0833 858 1 0.911 0.00929 0.893 0.929
## 7.1667 857 1 0.910 0.00934 0.892 0.928
## 7.2500 856 1 0.909 0.00939 0.890 0.927
## 7.3333 855 2 0.906 0.00949 0.888 0.925
## 7.5000 853 2 0.904 0.00959 0.886 0.923
## 7.5833 851 2 0.902 0.00968 0.883 0.921
## 7.7500 849 1 0.901 0.00973 0.882 0.920
## 7.8333 848 1 0.900 0.00978 0.881 0.919
## 8.0000 847 2 0.898 0.00987 0.879 0.918
## 8.0833 845 2 0.896 0.00996 0.877 0.916
## 8.1667 843 2 0.894 0.01005 0.874 0.914
## 8.3333 841 1 0.893 0.01009 0.873 0.913
## 8.4167 840 1 0.892 0.01013 0.872 0.912
## 8.5000 839 2 0.889 0.01022 0.870 0.910
## 8.5833 837 1 0.888 0.01026 0.869 0.909
## 8.6667 836 3 0.885 0.01039 0.865 0.906
## 8.7500 833 2 0.883 0.01047 0.863 0.904
## 8.8333 831 2 0.881 0.01056 0.861 0.902
## 8.9167 829 1 0.880 0.01060 0.859 0.901
## 9.0000 828 3 0.877 0.01072 0.856 0.898
## 9.0833 825 1 0.876 0.01076 0.855 0.897
## 9.3333 824 1 0.875 0.01080 0.854 0.896
## 9.4167 823 1 0.874 0.01083 0.853 0.895
## 9.5833 822 2 0.871 0.01091 0.850 0.893
## 9.6667 820 3 0.868 0.01103 0.847 0.890
## 9.7500 817 3 0.865 0.01114 0.843 0.887
## 9.8333 814 1 0.864 0.01118 0.842 0.886
## 10.0833 813 1 0.863 0.01121 0.841 0.885
## 10.2500 812 4 0.859 0.01136 0.837 0.881
## 10.5833 808 2 0.857 0.01143 0.834 0.879
## 10.6667 806 2 0.854 0.01150 0.832 0.877
## 10.7500 804 2 0.852 0.01157 0.830 0.875
## 10.9167 802 3 0.849 0.01167 0.827 0.872
## 11.0000 799 2 0.847 0.01174 0.824 0.870
## 11.0833 797 1 0.846 0.01177 0.823 0.869
## 11.1667 796 1 0.845 0.01180 0.822 0.868
## 11.2500 795 2 0.843 0.01187 0.820 0.866
## 11.3333 793 1 0.842 0.01190 0.819 0.865
## 11.4167 792 2 0.840 0.01197 0.816 0.863
## 11.5000 790 2 0.837 0.01203 0.814 0.861
## 11.5833 788 3 0.834 0.01212 0.811 0.858
## 11.6667 785 1 0.833 0.01215 0.810 0.857
## 11.7500 784 1 0.832 0.01218 0.809 0.856
## 11.8333 783 1 0.831 0.01222 0.807 0.855
## 11.9167 782 2 0.829 0.01228 0.805 0.853
## 12.0000 780 1 0.828 0.01231 0.804 0.852
## 12.0833 779 1 0.827 0.01234 0.803 0.851
## 12.1667 778 3 0.824 0.01243 0.800 0.848
## 12.2500 775 2 0.821 0.01248 0.797 0.846
## 12.3333 773 1 0.820 0.01251 0.796 0.845
## 12.4167 772 1 0.819 0.01254 0.795 0.844
## 12.5000 771 1 0.818 0.01257 0.794 0.843
## 12.6667 770 3 0.815 0.01266 0.791 0.840
## 12.7500 767 1 0.814 0.01268 0.790 0.839
## 12.8333 766 1 0.813 0.01271 0.788 0.838
## 13.0833 738 1 0.812 0.01274 0.787 0.837
## 13.1667 718 3 0.808 0.01284 0.784 0.834
## 13.2500 697 1 0.807 0.01287 0.782 0.833
## 13.3333 668 2 0.805 0.01295 0.780 0.831
## 13.4167 642 2 0.802 0.01303 0.777 0.828
## 13.5000 639 1 0.801 0.01307 0.776 0.827
## 13.5833 613 2 0.799 0.01315 0.773 0.825
## 13.6667 591 3 0.794 0.01329 0.769 0.821
## 13.8333 556 2 0.792 0.01340 0.766 0.818
## 14.0000 522 1 0.790 0.01346 0.764 0.817
## 14.0833 515 2 0.787 0.01358 0.761 0.814
## 14.2500 482 5 0.779 0.01392 0.752 0.807
## 14.3333 471 1 0.777 0.01399 0.750 0.805
## 14.4167 442 1 0.775 0.01407 0.748 0.804
## 14.6667 390 2 0.771 0.01427 0.744 0.800
## 15.0000 329 1 0.769 0.01442 0.741 0.798
## 15.2500 283 1 0.766 0.01463 0.738 0.796
## 15.6667 219 1 0.763 0.01497 0.734 0.793
## 15.7500 197 1 0.759 0.01539 0.729 0.790
## 15.9167 173 1 0.755 0.01591 0.724 0.786
## 16.1667 141 1 0.749 0.01668 0.717 0.783
## 16.5833 57 1 0.736 0.02093 0.696 0.778
## 16.6667 45 1 0.720 0.02609 0.670 0.773
##
## BOOZE_q=0.5–2/week
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0000 1729 1 0.999 0.000578 0.998 1.000
## 0.0833 1728 1 0.999 0.000817 0.997 1.000
## 0.1667 1727 1 0.998 0.001001 0.996 1.000
## 0.3333 1726 1 0.998 0.001155 0.995 1.000
## 0.4167 1725 1 0.997 0.001291 0.995 1.000
## 0.7500 1724 1 0.997 0.001414 0.994 0.999
## 0.8333 1723 1 0.996 0.001527 0.993 0.999
## 1.0000 1722 3 0.994 0.001824 0.991 0.998
## 1.1667 1719 2 0.993 0.001997 0.989 0.997
## 1.2500 1717 3 0.991 0.002230 0.987 0.996
## 1.3333 1714 2 0.990 0.002373 0.986 0.995
## 1.4167 1712 1 0.990 0.002441 0.985 0.994
## 1.5000 1711 1 0.989 0.002507 0.984 0.994
## 1.5833 1710 1 0.988 0.002572 0.983 0.993
## 1.6667 1709 2 0.987 0.002695 0.982 0.993
## 1.7500 1707 2 0.986 0.002814 0.981 0.992
## 1.8333 1705 1 0.986 0.002871 0.980 0.991
## 2.0000 1704 1 0.985 0.002927 0.979 0.991
## 2.0833 1703 2 0.984 0.003036 0.978 0.990
## 2.1667 1701 1 0.983 0.003088 0.977 0.989
## 2.2500 1700 1 0.983 0.003140 0.977 0.989
## 2.3333 1699 1 0.982 0.003191 0.976 0.988
## 2.5000 1698 4 0.980 0.003387 0.973 0.986
## 2.5833 1694 2 0.979 0.003480 0.972 0.985
## 2.6667 1692 2 0.977 0.003571 0.970 0.984
## 2.7500 1690 3 0.976 0.003702 0.968 0.983
## 2.8333 1687 3 0.974 0.003829 0.966 0.982
## 2.9167 1684 4 0.972 0.003991 0.964 0.980
## 3.0000 1680 1 0.971 0.004030 0.963 0.979
## 3.0833 1679 2 0.970 0.004107 0.962 0.978
## 3.1667 1677 2 0.969 0.004183 0.961 0.977
## 3.2500 1675 1 0.968 0.004221 0.960 0.976
## 3.3333 1674 2 0.967 0.004294 0.959 0.975
## 3.4167 1672 2 0.966 0.004366 0.957 0.974
## 3.5000 1670 1 0.965 0.004402 0.957 0.974
## 3.5833 1669 2 0.964 0.004472 0.955 0.973
## 3.6667 1667 3 0.962 0.004574 0.953 0.971
## 3.7500 1664 1 0.962 0.004608 0.953 0.971
## 3.8333 1663 3 0.960 0.004707 0.951 0.969
## 4.0000 1660 5 0.957 0.004868 0.948 0.967
## 4.0833 1655 3 0.955 0.004961 0.946 0.965
## 4.1667 1652 4 0.953 0.005082 0.943 0.963
## 4.2500 1648 1 0.953 0.005112 0.943 0.963
## 4.3333 1647 2 0.951 0.005170 0.941 0.962
## 4.4167 1645 1 0.951 0.005200 0.941 0.961
## 4.5000 1644 1 0.950 0.005228 0.940 0.961
## 4.5833 1643 1 0.950 0.005257 0.939 0.960
## 4.6667 1642 2 0.949 0.005314 0.938 0.959
## 4.7500 1640 2 0.947 0.005370 0.937 0.958
## 4.8333 1638 1 0.947 0.005398 0.936 0.957
## 4.9167 1637 2 0.946 0.005453 0.935 0.956
## 5.0000 1635 2 0.944 0.005507 0.934 0.955
## 5.0833 1633 4 0.942 0.005614 0.931 0.953
## 5.2500 1629 2 0.941 0.005666 0.930 0.952
## 5.4167 1627 2 0.940 0.005718 0.929 0.951
## 5.5000 1625 1 0.939 0.005744 0.928 0.951
## 5.5833 1624 1 0.939 0.005769 0.927 0.950
## 5.6667 1623 1 0.938 0.005795 0.927 0.950
## 5.7500 1622 2 0.937 0.005845 0.926 0.948
## 5.8333 1620 1 0.936 0.005870 0.925 0.948
## 6.0833 1619 2 0.935 0.005919 0.924 0.947
## 6.1667 1617 1 0.935 0.005944 0.923 0.946
## 6.2500 1616 1 0.934 0.005968 0.922 0.946
## 6.4167 1615 3 0.932 0.006041 0.921 0.944
## 6.5000 1612 1 0.932 0.006065 0.920 0.944
## 6.5833 1611 4 0.929 0.006159 0.917 0.942
## 6.6667 1607 3 0.928 0.006228 0.916 0.940
## 6.7500 1604 4 0.925 0.006319 0.913 0.938
## 6.8333 1600 1 0.925 0.006342 0.912 0.937
## 6.9167 1599 1 0.924 0.006364 0.912 0.937
## 7.0000 1598 1 0.924 0.006386 0.911 0.936
## 7.0833 1597 1 0.923 0.006408 0.911 0.936
## 7.1667 1596 1 0.922 0.006430 0.910 0.935
## 7.2500 1595 1 0.922 0.006452 0.909 0.935
## 7.3333 1594 4 0.920 0.006539 0.907 0.933
## 7.4167 1590 1 0.919 0.006560 0.906 0.932
## 7.5000 1589 1 0.918 0.006582 0.906 0.931
## 7.5833 1588 5 0.916 0.006687 0.903 0.929
## 7.6667 1583 6 0.912 0.006810 0.899 0.926
## 7.7500 1577 1 0.912 0.006830 0.898 0.925
## 7.8333 1576 1 0.911 0.006850 0.898 0.924
## 7.9167 1575 1 0.910 0.006870 0.897 0.924
## 8.0000 1574 3 0.909 0.006930 0.895 0.922
## 8.0833 1571 3 0.907 0.006989 0.893 0.921
## 8.1667 1568 4 0.905 0.007066 0.891 0.919
## 8.2500 1564 3 0.903 0.007123 0.889 0.917
## 8.3333 1561 3 0.901 0.007179 0.887 0.915
## 8.4167 1558 2 0.900 0.007217 0.886 0.914
## 8.5000 1556 3 0.898 0.007272 0.884 0.913
## 8.5833 1553 1 0.898 0.007290 0.883 0.912
## 8.6667 1552 1 0.897 0.007308 0.883 0.911
## 8.7500 1551 3 0.895 0.007363 0.881 0.910
## 8.8333 1548 1 0.895 0.007381 0.880 0.909
## 8.9167 1547 2 0.894 0.007416 0.879 0.908
## 9.0000 1545 1 0.893 0.007434 0.879 0.908
## 9.1667 1544 3 0.891 0.007487 0.877 0.906
## 9.4167 1541 1 0.891 0.007504 0.876 0.906
## 9.5000 1540 1 0.890 0.007521 0.875 0.905
## 9.5833 1539 3 0.888 0.007573 0.874 0.903
## 9.6667 1536 1 0.888 0.007590 0.873 0.903
## 9.7500 1535 2 0.887 0.007624 0.872 0.902
## 9.8333 1533 2 0.885 0.007658 0.871 0.901
## 9.9167 1531 2 0.884 0.007692 0.869 0.900
## 10.0000 1529 1 0.884 0.007708 0.869 0.899
## 10.0833 1528 3 0.882 0.007758 0.867 0.897
## 10.1667 1525 3 0.880 0.007807 0.865 0.896
## 10.2500 1522 2 0.879 0.007840 0.864 0.895
## 10.3333 1520 1 0.879 0.007856 0.863 0.894
## 10.4167 1519 4 0.876 0.007920 0.861 0.892
## 10.5000 1515 3 0.874 0.007967 0.859 0.890
## 10.5833 1512 5 0.872 0.008045 0.856 0.888
## 10.6667 1507 1 0.871 0.008061 0.855 0.887
## 10.7500 1506 1 0.870 0.008076 0.855 0.886
## 10.8333 1505 2 0.869 0.008107 0.854 0.885
## 10.9167 1503 2 0.868 0.008137 0.852 0.884
## 11.0000 1501 2 0.867 0.008167 0.851 0.883
## 11.1667 1499 1 0.866 0.008182 0.851 0.883
## 11.2500 1498 2 0.865 0.008212 0.849 0.881
## 11.3333 1496 2 0.864 0.008242 0.848 0.880
## 11.4167 1494 1 0.864 0.008256 0.847 0.880
## 11.5000 1493 3 0.862 0.008300 0.846 0.878
## 11.5833 1490 3 0.860 0.008344 0.844 0.877
## 11.6667 1487 2 0.859 0.008373 0.843 0.875
## 11.7500 1485 2 0.858 0.008401 0.841 0.874
## 11.8333 1483 3 0.856 0.008444 0.840 0.873
## 11.9167 1480 2 0.855 0.008472 0.838 0.872
## 12.0000 1478 2 0.854 0.008500 0.837 0.870
## 12.1667 1476 4 0.851 0.008555 0.835 0.868
## 12.2500 1472 3 0.850 0.008596 0.833 0.867
## 12.3333 1469 3 0.848 0.008637 0.831 0.865
## 12.4167 1466 3 0.846 0.008677 0.829 0.863
## 12.5000 1463 1 0.846 0.008690 0.829 0.863
## 12.6667 1462 2 0.844 0.008717 0.828 0.862
## 12.7500 1460 5 0.842 0.008782 0.824 0.859
## 12.8333 1455 2 0.840 0.008808 0.823 0.858
## 12.9167 1425 3 0.839 0.008849 0.821 0.856
## 13.0000 1403 1 0.838 0.008863 0.821 0.856
## 13.0833 1393 1 0.837 0.008877 0.820 0.855
## 13.1667 1364 4 0.835 0.008935 0.818 0.853
## 13.2500 1330 4 0.832 0.008996 0.815 0.850
## 13.3333 1287 2 0.831 0.009029 0.814 0.849
## 13.5000 1230 1 0.830 0.009046 0.813 0.848
## 13.5833 1182 1 0.830 0.009066 0.812 0.848
## 13.6667 1156 2 0.828 0.009107 0.811 0.846
## 13.7500 1119 2 0.827 0.009151 0.809 0.845
## 13.8333 1086 1 0.826 0.009174 0.808 0.844
## 13.9167 1071 1 0.825 0.009198 0.807 0.844
## 14.0833 1037 1 0.825 0.009223 0.807 0.843
## 14.1667 999 2 0.823 0.009278 0.805 0.841
## 14.2500 955 1 0.822 0.009308 0.804 0.840
## 14.3333 931 3 0.819 0.009403 0.801 0.838
## 14.4167 879 2 0.817 0.009474 0.799 0.836
## 14.5000 815 1 0.816 0.009515 0.798 0.835
## 14.9167 676 2 0.814 0.009639 0.795 0.833
## 15.0000 660 2 0.812 0.009766 0.793 0.831
## 15.0833 647 2 0.809 0.009896 0.790 0.829
## 15.1667 628 3 0.805 0.010097 0.786 0.825
## 15.5000 498 3 0.800 0.010418 0.780 0.821
## 15.5833 463 3 0.795 0.010772 0.774 0.817
## 15.9167 340 1 0.793 0.010991 0.772 0.815
##
## BOOZE_q=>2/week
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.167 2527 4 0.998 0.000791 0.997 1.000
## 0.250 2523 1 0.998 0.000884 0.996 1.000
## 0.333 2522 4 0.996 0.001185 0.994 0.999
## 0.417 2518 1 0.996 0.001249 0.994 0.998
## 0.500 2517 4 0.994 0.001477 0.992 0.997
## 0.583 2513 1 0.994 0.001528 0.991 0.997
## 0.750 2512 1 0.994 0.001578 0.991 0.997
## 0.833 2511 3 0.992 0.001718 0.989 0.996
## 0.917 2508 1 0.992 0.001763 0.989 0.996
## 1.000 2507 2 0.991 0.001848 0.988 0.995
## 1.083 2505 4 0.990 0.002007 0.986 0.994
## 1.167 2501 2 0.989 0.002082 0.985 0.993
## 1.250 2499 1 0.989 0.002119 0.984 0.993
## 1.333 2498 3 0.987 0.002224 0.983 0.992
## 1.417 2495 1 0.987 0.002258 0.983 0.991
## 1.500 2494 3 0.986 0.002357 0.981 0.990
## 1.583 2491 1 0.985 0.002389 0.981 0.990
## 1.750 2490 1 0.985 0.002421 0.980 0.990
## 1.833 2489 2 0.984 0.002483 0.979 0.989
## 2.000 2487 1 0.984 0.002513 0.979 0.989
## 2.083 2486 1 0.983 0.002543 0.978 0.988
## 2.167 2485 1 0.983 0.002573 0.978 0.988
## 2.250 2484 1 0.983 0.002602 0.978 0.988
## 2.333 2483 1 0.982 0.002631 0.977 0.987
## 2.417 2482 1 0.982 0.002659 0.977 0.987
## 2.500 2481 4 0.980 0.002770 0.975 0.986
## 2.583 2477 2 0.979 0.002824 0.974 0.985
## 2.667 2475 1 0.979 0.002851 0.973 0.985
## 2.833 2474 3 0.978 0.002928 0.972 0.984
## 2.917 2471 3 0.977 0.003004 0.971 0.983
## 3.000 2468 6 0.974 0.003149 0.968 0.980
## 3.083 2462 3 0.973 0.003219 0.967 0.979
## 3.167 2459 1 0.973 0.003242 0.966 0.979
## 3.250 2458 3 0.972 0.003310 0.965 0.978
## 3.333 2455 2 0.971 0.003354 0.964 0.977
## 3.417 2453 1 0.970 0.003376 0.964 0.977
## 3.500 2452 3 0.969 0.003441 0.962 0.976
## 3.583 2449 3 0.968 0.003504 0.961 0.975
## 3.667 2446 2 0.967 0.003546 0.960 0.974
## 3.750 2444 1 0.967 0.003566 0.960 0.974
## 3.833 2443 3 0.966 0.003627 0.958 0.973
## 3.917 2440 2 0.965 0.003667 0.958 0.972
## 4.000 2438 4 0.963 0.003745 0.956 0.971
## 4.250 2434 1 0.963 0.003765 0.955 0.970
## 4.333 2433 1 0.962 0.003784 0.955 0.970
## 4.417 2432 2 0.962 0.003822 0.954 0.969
## 4.500 2430 3 0.960 0.003878 0.953 0.968
## 4.583 2427 2 0.960 0.003915 0.952 0.967
## 4.667 2425 4 0.958 0.003988 0.950 0.966
## 4.750 2421 4 0.956 0.004059 0.949 0.964
## 4.833 2417 3 0.955 0.004111 0.947 0.963
## 4.917 2414 2 0.954 0.004146 0.946 0.963
## 5.000 2412 4 0.953 0.004214 0.945 0.961
## 5.083 2408 6 0.951 0.004314 0.942 0.959
## 5.167 2402 1 0.950 0.004330 0.942 0.959
## 5.250 2401 4 0.949 0.004394 0.940 0.957
## 5.333 2397 1 0.948 0.004410 0.940 0.957
## 5.417 2396 4 0.947 0.004473 0.938 0.955
## 5.500 2392 2 0.946 0.004505 0.937 0.955
## 5.583 2390 3 0.945 0.004551 0.936 0.954
## 5.667 2387 3 0.943 0.004596 0.934 0.952
## 5.750 2384 7 0.941 0.004701 0.931 0.950
## 5.833 2377 5 0.939 0.004773 0.929 0.948
## 5.917 2372 4 0.937 0.004830 0.928 0.947
## 6.000 2368 1 0.937 0.004845 0.927 0.946
## 6.167 2367 2 0.936 0.004873 0.926 0.945
## 6.250 2365 2 0.935 0.004901 0.926 0.945
## 6.333 2363 5 0.933 0.004969 0.923 0.943
## 6.417 2358 2 0.932 0.004997 0.923 0.942
## 6.583 2356 3 0.931 0.005037 0.921 0.941
## 6.667 2353 8 0.928 0.005143 0.918 0.938
## 6.750 2345 2 0.927 0.005169 0.917 0.937
## 6.833 2343 1 0.927 0.005182 0.917 0.937
## 6.917 2342 1 0.926 0.005195 0.916 0.937
## 7.000 2341 4 0.925 0.005246 0.915 0.935
## 7.083 2337 2 0.924 0.005271 0.914 0.934
## 7.167 2335 5 0.922 0.005333 0.912 0.933
## 7.250 2330 4 0.920 0.005383 0.910 0.931
## 7.333 2326 4 0.919 0.005431 0.908 0.930
## 7.417 2322 3 0.918 0.005467 0.907 0.928
## 7.500 2319 3 0.917 0.005503 0.906 0.927
## 7.583 2316 7 0.914 0.005585 0.903 0.925
## 7.667 2309 7 0.911 0.005665 0.900 0.922
## 7.750 2302 2 0.910 0.005688 0.899 0.921
## 7.833 2300 3 0.909 0.005722 0.898 0.920
## 8.000 2297 3 0.908 0.005755 0.897 0.919
## 8.083 2294 3 0.907 0.005788 0.895 0.918
## 8.167 2291 2 0.906 0.005810 0.895 0.917
## 8.250 2289 1 0.905 0.005821 0.894 0.917
## 8.333 2288 4 0.904 0.005865 0.892 0.915
## 8.417 2284 4 0.902 0.005908 0.891 0.914
## 8.500 2280 7 0.899 0.005981 0.888 0.911
## 8.583 2273 2 0.899 0.006002 0.887 0.911
## 8.667 2271 5 0.897 0.006054 0.885 0.909
## 8.750 2266 4 0.895 0.006095 0.883 0.907
## 8.833 2262 6 0.893 0.006155 0.881 0.905
## 8.917 2256 5 0.891 0.006205 0.879 0.903
## 9.000 2251 5 0.889 0.006254 0.877 0.901
## 9.083 2246 4 0.887 0.006293 0.875 0.900
## 9.167 2242 4 0.886 0.006331 0.873 0.898
## 9.250 2238 2 0.885 0.006350 0.872 0.897
## 9.333 2236 2 0.884 0.006369 0.872 0.897
## 9.417 2234 4 0.882 0.006407 0.870 0.895
## 9.500 2230 4 0.881 0.006444 0.868 0.894
## 9.583 2226 6 0.879 0.006499 0.866 0.891
## 9.667 2220 3 0.877 0.006526 0.865 0.890
## 9.750 2217 3 0.876 0.006553 0.863 0.889
## 9.833 2214 3 0.875 0.006580 0.862 0.888
## 9.917 2211 1 0.875 0.006589 0.862 0.888
## 10.000 2210 2 0.874 0.006607 0.861 0.887
## 10.083 2208 5 0.872 0.006651 0.859 0.885
## 10.167 2203 5 0.870 0.006694 0.857 0.883
## 10.250 2198 2 0.869 0.006712 0.856 0.882
## 10.333 2196 3 0.868 0.006737 0.855 0.881
## 10.417 2193 7 0.865 0.006797 0.852 0.878
## 10.500 2186 3 0.864 0.006822 0.851 0.877
## 10.583 2183 1 0.863 0.006830 0.850 0.877
## 10.667 2182 6 0.861 0.006880 0.848 0.875
## 10.750 2176 3 0.860 0.006904 0.846 0.874
## 10.833 2173 2 0.859 0.006921 0.846 0.873
## 10.917 2171 5 0.857 0.006961 0.844 0.871
## 11.000 2166 6 0.855 0.007009 0.841 0.869
## 11.083 2160 4 0.853 0.007040 0.839 0.867
## 11.167 2156 1 0.853 0.007048 0.839 0.867
## 11.250 2155 3 0.852 0.007072 0.838 0.866
## 11.333 2152 2 0.851 0.007087 0.837 0.865
## 11.417 2150 3 0.850 0.007110 0.836 0.864
## 11.500 2147 3 0.848 0.007134 0.835 0.863
## 11.583 2144 3 0.847 0.007156 0.833 0.861
## 11.667 2141 6 0.845 0.007202 0.831 0.859
## 11.750 2135 4 0.843 0.007232 0.829 0.858
## 11.833 2131 2 0.843 0.007246 0.828 0.857
## 11.917 2129 1 0.842 0.007254 0.828 0.856
## 12.000 2128 1 0.842 0.007261 0.828 0.856
## 12.083 2127 2 0.841 0.007276 0.827 0.855
## 12.167 2125 5 0.839 0.007312 0.825 0.853
## 12.250 2120 4 0.837 0.007341 0.823 0.852
## 12.333 2116 3 0.836 0.007363 0.822 0.851
## 12.417 2113 8 0.833 0.007419 0.819 0.848
## 12.500 2105 3 0.832 0.007441 0.817 0.847
## 12.583 2102 2 0.831 0.007454 0.817 0.846
## 12.667 2100 5 0.829 0.007489 0.814 0.844
## 12.750 2095 1 0.829 0.007496 0.814 0.843
## 12.833 2094 4 0.827 0.007523 0.812 0.842
## 12.917 2060 2 0.826 0.007537 0.812 0.841
## 13.000 2045 3 0.825 0.007559 0.810 0.840
## 13.083 2032 4 0.823 0.007587 0.809 0.838
## 13.167 1988 5 0.821 0.007625 0.807 0.836
## 13.250 1948 2 0.821 0.007640 0.806 0.836
## 13.333 1886 6 0.818 0.007690 0.803 0.833
## 13.417 1804 7 0.815 0.007753 0.800 0.830
## 13.500 1794 3 0.813 0.007780 0.798 0.829
## 13.583 1720 2 0.812 0.007799 0.797 0.828
## 13.667 1681 5 0.810 0.007851 0.795 0.826
## 13.750 1646 4 0.808 0.007893 0.793 0.824
## 13.833 1612 4 0.806 0.007937 0.791 0.822
## 13.917 1589 3 0.805 0.007970 0.789 0.820
## 14.000 1565 1 0.804 0.007982 0.789 0.820
## 14.083 1538 1 0.803 0.007994 0.788 0.819
## 14.167 1464 2 0.802 0.008020 0.787 0.818
## 14.250 1414 3 0.801 0.008063 0.785 0.817
## 14.333 1387 1 0.800 0.008078 0.784 0.816
## 14.417 1328 2 0.799 0.008111 0.783 0.815
## 14.500 1260 5 0.796 0.008202 0.780 0.812
## 14.583 1220 2 0.794 0.008240 0.778 0.811
## 14.667 1188 1 0.794 0.008260 0.778 0.810
## 14.750 1115 4 0.791 0.008352 0.775 0.807
## 14.833 1041 1 0.790 0.008379 0.774 0.807
## 14.917 1001 1 0.789 0.008407 0.773 0.806
## 15.000 963 2 0.788 0.008469 0.771 0.804
## 15.333 808 2 0.786 0.008560 0.769 0.803
## 15.417 755 3 0.783 0.008714 0.766 0.800
## 15.583 614 3 0.779 0.008946 0.761 0.797
## 15.667 589 1 0.777 0.009028 0.760 0.795
## 15.750 548 1 0.776 0.009123 0.758 0.794
## 15.833 530 2 0.773 0.009320 0.755 0.792
## 15.917 492 1 0.772 0.009433 0.753 0.790
## 16.000 460 1 0.770 0.009560 0.751 0.789
## 16.083 441 1 0.768 0.009697 0.749 0.787
## 16.167 422 1 0.766 0.009843 0.747 0.786
## 16.250 362 3 0.760 0.010422 0.740 0.781
## 16.417 225 1 0.757 0.010910 0.736 0.778
## 16.500 187 2 0.749 0.012201 0.725 0.773
## 16.583 169 1 0.744 0.012908 0.719 0.770
## 16.750 82 1 0.735 0.015618 0.705 0.766
summary(fit)$table
## records n.max n.start events rmean se(rmean) median
## BOOZE_q=0/week 4053 4053 4053 1070 14.42166 0.07071389 NA
## BOOZE_q=0–0.5/week 941 941 941 213 14.83118 0.13643836 NA
## BOOZE_q=0.5–2/week 1729 1729 1729 325 15.15602 0.09331554 NA
## BOOZE_q=>2/week 2527 2527 2527 537 15.03836 0.07729027 NA
## 0.95LCL 0.95UCL
## BOOZE_q=0/week NA NA
## BOOZE_q=0–0.5/week NA NA
## BOOZE_q=0.5–2/week NA NA
## BOOZE_q=>2/week NA NA
#Log-rank test
survdiff(Surv(FU, DEATH)~BOOZE_q, data=d)
## Call:
## survdiff(formula = Surv(FU, DEATH) ~ BOOZE_q, data = d)
##
## N Observed Expected (O-E)^2/E (O-E)^2/V
## BOOZE_q=0/week 4053 1070 916 25.873 45.243
## BOOZE_q=0–0.5/week 941 213 219 0.139 0.155
## BOOZE_q=0.5–2/week 1729 325 411 18.013 22.319
## BOOZE_q=>2/week 2527 537 599 6.493 9.027
##
## Chisq= 50.6 on 3 degrees of freedom, p= 6e-11